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Abstract 

We compute perturbative QCD corrections to the lifetime splitting between the charged and 
neutral B meson in the framework of the heavy quark expansion. These next-to-leading log- 
arithmic corrections are necessary for a meaningful use of hadronic matrix elements of lo- 
cal operators from lattice gauge theory. We find the uncertainties associated with the choices 
of renormalization scale and scheme significantly reduced compared to the leading-order re- 
sult. We include the full dependence on the charm-quark mass m c without any approxima- 
tions. Using hadronic matrix elements estimated in the literature with lattice QCD we obtain 
t(B + )/t(B®) = 1.053 ± 0.016 ± 0.017, where the effects of unquenching and l/m b correc- 
tions are not yet included. The lifetime difference of heavy baryons S° and 2^ is also briefly 
discussed. 
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Figure 1: Weak annihilation (WA) and Pauli interference (PI) diagrams in the leading order 
of QCD. They contribute to T(B^) and T(B + ), respectively. The crosses represent \AB\ = 1 
operators, which are generated by the exchange of W bosons. CKM-suppressed contributions 
are not shown. 

1 Preliminaries 

The Heavy Quark Expansion (HQE) technique provides a well-defined QCD-based framework 
for the calculation of total decay rates of b- flavoured hadrons [[[]]. The HQE yields an expansion 
of the decay rate T(H b ) in terms of AQ C£) /m b , where H b represents any hadron containing a 
single 6-quark and any of the light u,d,s (anti-)quarks as valence quarks. m b is the 6-quark mass 
and Aqcd is the fundamental scale of QCD, which determines the size of hadronic effects. In 
the leading order of h.QCD/ m b the decay rate of H b equals the decay rate of a free 6-quark, 
which is unaffected by the light degrees of freedom of H b . Consequently, the lifetimes of all 
6-flavoured hadrons are the same at this order. The first corrections to the free quark decay 
appear at order (Ag^/m^) 2 an d are caused by the Fermi motion of the 6-quark in H b and the 
chromomagnetic interaction of the final state quarks with the hadronic cloud surrounding the 
heavy 6-quark. These mechanisms have a negligible effect on the lifetime difference between 
the B + and B® mesons, because the strong interaction excellently respects isospin symmetry. At 
order (Ag^/mf,) 3 , however, one encounters weak interaction effects between the 6-quark and 
the light valence quark. These effects, known as weak annihilation (WA) and Pauli interference 
(PI) [Jjl, are depicted in Fig. [I]. They are phase-space enhanced with respect to the leading 
free-quark decay and induce corrections to T(H b ) of order 167r 2 (AgcD/mf,) 3 = 0(5—10%). 
The measurement of lifetime differences among different 6-flavoured hadrons therefore tests the 
HQE formalism at the third order in the expansion parameter. 

The calculation of T(H b ) consists of three steps: the first step is an operator product expan- 
sion (OPE) integrating out the heavy W boson, which mediates the weak b decay. This results in 
an effective \AB\ = 1 Hamiltonian describing the flavour-changing weak interaction of the Stan- 
dard Model up to corrections of order rr? b jMy V , where AB denotes the change in bottom-quark 
number: 



G 



H = E v wd> Cx^OQf^i) + C 2 0i)Q|' d '0i) + h.c. (1) 



V2 c \, 



=d,s 



Here Gp is the Fermi constant and the V^ 's are elements of the Cabibbo-Kobayashi-Maskawa 
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(CKM) matrix. The Wilson coefficients C,(//i) contain the short-distance physics associated 
with scales above the renormalization scale The weak interaction is encoded in the four- 
quark operators 



Q\ d =^(i- 75 ) Cj ^(i- 75 )^, 



6i 7M (l- 75)^^(1 -75)4, (2) 



with summation over the colour indices i and j. We have omitted penguin operators and doubly 
Cabibbo-suppressed terms in (P, which have a negligible effect on the B + -B® lifetime differ- 
ence. Next the total decay rate T(H b ) is related to H by the optical theorem: 



T(H b 



2M 



■(H b \T\H b ). 



(3) 



Here we have adopted the conventional relativistic normalization (H b \H b ) = 2EV and intro- 
duced the transition operator: 



T 



Imi drxT\H(x)H 



(4) 



The second step is the HQE, which exploits the hierarchy m b ^> Aqcd to expand the RHS of (Q) 
in terms of Aqct)/ 77 ^- To this end an OPE is applied to T which effectively integrates out the 
hard loop momenta (corresponding to the momenta of the final state quarks). We decompose the 
result as 



r 



[T + T 2 + r 3 
T u + T d + T s 



1 



0(l/m 



sing 



(5) 



Here T n denotes the portion of T which is suppressed by a factor of l/m b with respect to % 
describing the free quark decay. The contributions to % from weak spectator interactions read 



T 



G 2 F m 2 b \V cb \ 2 

67T 



G 2 F m 2 b \V cb \ 2 
6ti 



\V ud \ 2 (F u Q d + FgQ d s + G u T d + G U S T C S 
+ \V cd \ 2 [F c Q d + F c s Q d s + G c T d + G C S T, 



+ {d 



F d Q u + F£Q U S + G a T u + G a s Tg 



(6) 



The superscript q of the coefficients F q , Ff, G q , G q s refers to the cq intermediate state (see 
Fig. []]). We include singly Cabibbo-suppressed contributions. In writing T d we have used 
\Vud\ 2 + |^4s| 2 ~ 1 an d m d ~ m s ~ 0, so that F d = F s , etc.. Here we encounter the local 
dimension-6, AB = operators 



Q q = b-fJl-^QQ^d 



T q = h»(l-l 5 )T a qqY(l-l 5 )T a b, 



Q q 



1 S 



s - &( 1_ 75)??(l+75)6, 

b(l- lFj )T a qq(l + l5 )T a b, 



(7) 
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where T a is the generator of colour SU(3). We define the AB = operators at the renormaliza- 
tion scale /i , which is of order m b . The Wilson coefficients F u . . . G d s are computed in pertur- 
bation theory. When applied to mesons, T u and T d correspond to the WA and PI mechanisms 
of Fig. [TJ respectively. In the case of baryons their role is interchanged: T u encodes the PI effect 
and T d describes the weak scattering of the 6-quark with the valence quark (see Fig. |5|). The co- 
efficients in @ depend on fi . Since the hard loops involve the charm quark, they also depend on 
the ratio z = m 2 /ml. The truncation of the perturbation series makes F u . . . G d s also dependent 
on \i\ = 0(rtib). This dependence diminishes in increasing orders of a s . To the considered or- 
der, the dependence on /z cancels between the coefficients and the matrix elements of operators 
in (0), so that observables are independent of /i . The remainder T sing in (|5|) involves additional 
dimension-6 operators, which describe power-suppressed contributions to the free quark decay 
from strong interactions with the spectator quark. The operators in T sing are isospin singlets and 
do not contribute to the B + -B® lifetime difference. The formalism of (@)-(0) applies to weakly 
decaying hadrons containing a single bottom quark and no charm quarks. Decays of hadrons like 
the B c meson with more than one heavy quark have a different power counting than in (§) ||2p. In 
the third step one computes the hadronic matrix elements of the operators in ([7]). They enter our 
calculation in isospin-breaking combinations and are conventionally parametrized as [pj] 



(B + \(Q U 
{B + \(T U 



Q d )(Vo)\B + ) 



f B MlB x {^) 
f B M 2 B eM, 



(B + \(Q u s -Q d s)(Vo)\B + ) 

(B + m-Ti)(» )\B+) ■. 



: flMlBM, 

f 2 B M 2 B eM- (8) 



Here f B is the B meson decay constant. In the vacuum saturation approximation (VSA) one has 
i?i(/i ) = 1, B 2 {fi ) = 1 + 0(a s (m b ),A Q cD/m b ) and 61,2(^0) = 0. Corrections to the VSA 
results are of order 1 /N c , where N c = 3 is the number of colours. 

Using the isospin relation (B° d \Q d ' u \B%) = (B + \Q u ' d \B + ) we now find from © and (§): 



r(i?2) - r(B" 



12tt 



ud I 



2 F U + \V cd \ 2 F c -F d ) ■ B. 



(9) 



Here we have introduced the shorthand notation 



( F*(z,fM>) \ 
F q s (z,fio) 
G«(z,no) 



B(nc 



( B 1 M \ 

B 2 (no) 
e i(/^o) 



for q = d,u, c. (10) 



The strong interaction affects all three steps of the calculation. The minimal way to in- 
clude QCD effects is the leading logarithmic approximation, which includes corrections of order 
a^\n n {p,i/M w ), n — 0, 1, ... in the coefficients Ci )2 (a«i) in (jl]). The corresponding leading or- 
der (LO) calculation of the width difference in (Q) involves the diagrams in Fig. |T] [[lJBJ]. Yet LO 
results are too crude for a precise calculation of lifetime differences. The heavy-quark masses in 
(Q) cannot be defined in a proper way and one faces a large dependence on the renormalization 
scale fjLi. Furthermore, results for £?i )2 and e\ j2 from lattice gauge theory cannot be matched to the 
continuum theory in a meaningful way at LO. Finally, as pointed out in [[3j], at LO the coefficients 
F, Fg in <Q) are anomalously small. They multiply the large matrix elements parametrized by 
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£>i j2 , while the larger coefficients G, Gs come with the small hadronic parameters ei 2 , rendering 
the LO prediction highly unstable. To cure these problems one must include the next-to-leading- 
order (NLO) QCD corrections of order \n n (ni/M w ). NLO corrections to the effective 
\AB\ = 1 Hamiltonian in (|l]) have been computed in The second step beyond the LO re- 

quires the calculation of QCD corrections to the coefficients F u . . . Gf, in (Q). Such a calculation 
has been first performed for the B®-B® lifetime difference in where 0(a s ) corrections were 
calculated in the SU(3)f limit neglecting certain terms of order z. In this limit only a few penguin 
effects play a role. A complete NLO computation has been carried out for the lifetime difference 
between the two mass eigenstates of the B® meson in [|7|]. In particular the correct treatment 
of infrared effects, which appear at intermediate steps of the calculation, has been worked out 
in [^]]. The computation presented in this paper is conceptually similar to the one in [f7|], except 
that the considered transition is AB = rather than AB = 2 and the quark masses in the final 
state are different. While this work was in preparation, QCD corrections to T u and T d have also 
been calculated in [H], There are two important differences between our analysis and [|8p: 



(i) in [||] the NLO corrections have been computed for the limiting case z = 0, i.e. neglecting 
the charm-quark mass in the final state. The corrections to this limit are of order z In z or 
roughly 20%. In Sect. [2] we include the dependence on the charm-quark mass exactly. 



(ii) in [§| the AB = operators have been defined in the heavy quark effective theory (HQET) 
rather than in full QCD, as we did in ([7]). HQET operators were chosen to eliminate 
the mixing of the dimension-6 operators in (0) into lower-dimensional operators under 
renormalization. We emphasize that this mixing does not impede the use of QCD operators 
in the HQE: it results purely from ultraviolet effects and can be accounted for by a finite 
renormalization of the affected operators. For a more detailed discussion with an explicit 



example we refer the reader to P] and to Sect. 3.2 



Finally one must compute the non-perturbative QCD effects residing in f%Bi, . . . /^ e 2- Results 
from lattice gauge theory for the matrix elements in have been recently obtained in [Q]. Earlier 
results using HQET fields can be found in []T0[]. In the matching of the results to continuum QCD 
the dependence of B%, ... e-i on fi and on the chosen renormalization scheme must cancel the 
corresponding dependence of the Wilson coefficients, which requires NLO accuracy. 

2 T u and T d at next-to-leading order 

We decompose the Wilson coefficients in @ as 

F u {z,fi ) = CfMF&fax^Xn) + C 1 {fi 1 )C 2 {fii)F^{z,x fll ,x^) 

F^x^xJ = F<°\z) +^^F< 1 \z,x, 1 ,x^) + 0(al) (11) 
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with x M = fi/rrib and an analogous notation for the remaining Wilson coefficients in @. The LO 
coefficients are obtained from the diagrams in Fig. [I]. The non-vanishing coefficients read 0] 

= - (1 - zf (l + | 



1 pu,(0) 

3 11 


(*) 


_ 1 p«,(0) 
- ^12 


(*) 


= 3J#°>(z) 


1 Wo) 

3*5,11 
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- 2^5,12 
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= 3F 5 u j?(z) 
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= i*fi (0) ( 
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= 3F 2 c 2 (0) (z): 






1 P c,(0), 
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= 3*g£ ) (*)- 


6F d,(0) ( 






) = 


6F 2 d 2 (0) (z) = 



= -G^(^) = (i-zr (i+2,), 

^ 2 (0) (z) = - v / T^4i(l-z), (12) 



1 



-^G c ^>(z) = VT^(l + 2z), 
Gt?\z)=Gtt \z)=6(l-zf, 



while 



^,(0) _ ^,(0) _ ^,(0) _ ^,(0) _ ^0,(0) _ n c,{Q) _ n c,(0) _ n c,(0) _ n d,(0) _ n 
^11 — ^12 — ^5,11 ~~ ^5,12 — ^11 — ^12 — ^5,11 — Lt 5,12 ~~ ^12 — u ; 

F*f = G s f = 0. (13) 

To obtain the NLO corrections F^'^ . . . G s '^ we have calculated the diagrams Ei and the 
imaginary parts of Di in Fig. [Z[ At NLO one becomes sensitive to the renormalization scheme. 
First, this affects the quantities m&, z and a s entering our calculation. The NLO coefficients 
given below correspond to the use of the pole-mass definition for m& and the definition of a s in 
the MS scheme [|TT|]. z can be either calculated from the pole masses or from the MS masses, 
because z = m 2 c /m 2 = ml(m c ) /ml(m b ) + 0{a 2 s ). Second, the choice of the renormalization 
scheme is also an issue for the effective four-quark operators appearing at the various stages of 
our calculation. In the prediction of physical quantities this scheme dependence cancels to the 
calculated order, nevertheless it must be taken care of when assembling pieces from different 
theoretical sources. The Wilson coefficients C 1)2 of H in ([!]) and F™'^ . . . Gg^ depend on 
the scheme used to renormalize the AB = 1 operators in (Q), but this dependence cancels in 
F M 'W . . . Gg^\ Our results below correspond to the definition of Ci )2 in [^]. F u '^ . . . Gg'^ 1 ' 
also depend on the renormalization scheme of the AB = operators in (0). This dependence 
cancels only when these coefficients are combined with the hadronic parameters £>i 2 and 
calculated from lattice QCD. It is therefore important that our scheme is used in the lattice- 
continuum matching of these quantities. We use the MS scheme with the NDR prescription for 
75 ©• To specify the scheme completely, it is further necessary to state the definition of the 
evanescent operators appearing in the calculation [|T2|]. We use 

E[Q] = h, lplu (l- l5 )qqYl P r(l-l 5 )b - (4 - 8e) Q 
E[Q S ] = 67^(1 -75)??7V(1 + 75)6 - (4-8e)Q s (14) 



and analogous definitions of E[T] and E[T S }. When the diagrams Ei . . . E4 for e.g. Q$ are 
calculated in D = 4 — 2e dimensions, the result can be expressed as a linear combination of 
Qs and E[Qs). Effectively, (jT3|) defines how Dirac strings with two or three Dirac matrices are 
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Figure 2: WA contributions in the next-to-leading order of QCD. The PI diagrams are obtained 
by interchanging u and d and reversing the fermion flow of the u and d lines. The first line 
shows the radiative corrections to AB = operators, which are necessary for the proper infrared 
factorization. Not displayed are the diagrams E' 3 , E4 and D 3 _ 8 which are obtained from the 
corresponding unprimed diagrams by left-right reflection and the reverse of the fermion flow. 



reduced. (Note that < JT4| ) also implies the replacement rules 7a,7 p 7 At (l — 75) <E> 7 !y 7 p 7 M (l — 75) — ► 
(16-4£)7^(l-7 5 )®7^(l-7 5 )and7^(l-7 5 )®7'V(l+7 5 ) -> 4(l + £)(l- 75 )®(l + 7 5 ).) 
The particular choice of the — 8e terms in ([14]) is motivated by Fierz invariance: the one-loop 
matrix elements of e.g. Q s and its Fierz transform = —1/2 6^7^(1 + 7 5 )&jg 7 -7 ly (l — 75)^ 
are in general different. This feature is an artifact of dimensional regularization. With ( fT4| ) and 
a corresponding definition of E[Qg], however, Fierz invariance is maintained at the one-loop 
level. This choice, which has also been made in [Q] for the AB = 1 operators, has the practical 
advantage that one can freely use the Fierz transformation at any step of the calculation. In other 
words: "Fierz-evanescent" operators like Qs — Q$ can be identified with 0. 

In the procedure of matching the full theory (eq. (Q)) to the effective AB = theory, infrared 
singularities are encountered at 0(a s ) both in the full-theory diagrams and in the matrix elements 
of operators in the effective theory. The diagrams relevant for this issue are Di - _D 4 and E 1 - 
£4. The singularities cancel in the Wilson coefficients F and G, but need to be regularized at in- 
termediate steps of the calculation. We take the 6-quark on-shell, assign zero 4-momentum to the 
external light quarks and use dimensional regularization for the infrared (as well as the ultravio- 
let) divergences. In this case, care has to be taken to treat the Dirac algebra in a consistent way. 
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In computing the matching condition between D\ - D4 and E% - E4 we have used two different 
methods, which lead to the same result. In both methods ultraviolet divergences appearing in E x 
- E4 and D 3 are subtracted, respectively, by AB = and AB = 1 counterterms, in the usual 
way. 

In the first method, we distinguish IR singularities arising in loop integrals from UV sin- 
gularities, and treat the Dirac algebra in strictly four dimensions in the IR-divergent parts. In 
the second method, IR and UV divergences are not distinguished and (/-dimensional Dirac al- 
gebra is used throughout. In this case evanescent operators E, as those given in (|14|), give a 
non-vanishing contribution in the matching procedure. This is a subtlety of the IR regulator used 
in method 2 [O]. If a different IR regulator, such as a gluon mass or method 1, is used, the 
non-vanishing bare one-loop matrix element of E is cancelled by a finite counterterm, so that E 
disappears from the NLO matching calculation fl5|,|I2p. The non-zero contribution in method 2 
originates in diagram E\ with the insertion of an evanescent operator E. This diagram is zero 
in dimensional regularization, thus leaving the corresponding counterterm uncancelled. We have 
further parametrized the evanescent 0(e) parts appearing in the d-dimensional projections of 
general Dirac structures r (g> T onto the basic operators Q and Qs- There are four independent 
parameters in the calculation, corresponding to T being a string of two, three, four or five Dirac 
matrices. We have checked that all four parameters disappear from the final result for the coeffi- 
cients. (This is true for the evanescent 0(e) parts multiplying IR poles. The UV poles give rise 
to a dependence on these parameters, which corresponds to a usual scheme dependence that is 
cancelled by the matrix elements of operators in the effective theory. Our choice of scheme is 
specified by ([14]).) 

We would also like to mention that the Fierz ordering of AB = 1 operators is immaterial 
because Fierz symmetry is respected by the standard NDR renormalization scheme employed by 
us. This has been checked by using the Fierz form leading to Dirac strings with flavour structure 
bb ® uu in method 1, and bu ® ub in method 2, and similarly for the contribution with u — > d. 
(The Fierz form used in method 2 for bd ® db is such that a closed fermion loop is generated in 
Di - D 4 .) 

In the NLO corrections to (Q) we set \ V u d\ = 1 and V C( i = 0. This introduces an error of order 
\Vcd\ 2 a s (mb)zlnz, which is well below 1% of t(B + )/t(B^) — 1. Hence (g) only involves the 
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2 (1 -z) (60 + 77 z + 7 z 2 ) 
9 



2 (1 - z) (-2803 + 2786 z + 725 z 2 ) 
71 + 27 



^s'n ( z i x lJ-i ' X W ) ^5,11 ( z i x m 1 x IJ-o ) ~~ 



-18 (l-*) 2 (l + 2z)] ln(l-z) + 

4 (1-z) (28 + 103z-164^ 2 ) 
9 



-44 (4 - 3 z) z 1 



Mz) + 



Gs'fi2 i z ) x Hi i X IJ-o ) ^s}l2 i z i x Vi ' X Vo ) 



[16 (1-z) (l + 2z)} 



U 2 (z) + 



ln(l - z) \n(z) 



+ 



-4(1-^) 2 (1 + z) (l + 2z) 



ln(l-z) + 



Az (6-51z + 28z 2 



\n(z) + 
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-24 (1-z) 2 {l + 2z)\ ln(z w ) + 

4 (I-2) (-53 - 80z±82z 2 ) 
9 



[1-z) z (l±2z) 



7T 2 + 



16 (1-z) (l±2z) (3±5z " 



Li 2 (z) 



-2(l-zr (-2 + 31z±64z 2 ±3z 3 



3z 

2 z (36 - 336 2 + 62 z 2 + 9 z 3 



ln(l - z) ln(z) 



ln(l-z) + 



+ 



9 



[1-zY (l + 2z) 



ln(a; w ) + 



ln(z) + 



4(l-zr (l + 2z) 



ln(x M0 ) + 



-4 (1-z) (l+_2z) (9±7z) 
9 



2 (1-z) (385 + 1519 z - 3278 z 2 

7T + 

27 



(15) 



Here Li 2 (z) = — J 2 dt [ln(l — £)]/£ is the dilogarithm function. Any dependence on infrared 
regulators has cancelled from the coefficients in ( |T3| ) showing that infrared effects properly fac- 
torize. As another check we have verified that the dependence on ix\ cancels analytically to the 
calculated order. 

For our numerical studies we choose the following range for the input parameters: 



a s {M z 



0.118 ±0.003, 



4.8±0.1GeV, 



0.085 ±0.015. (16) 



Throughout this paper we always remove 0(a 2 s ) terms from the calculated coefficients. (For 
instance, at NLO we write a product such as C\F U as Cilo-^nlo + 2Ci,lo dC\F^ , where 
Ci.nlo = Ci,lo + dCi denotes the NLO Wilson coefficient.) In all terms we use the two-loop 
expression for the running coupling a s in QCD with five flavours. Numerical values for the 
calculated coefficients can be found in Tab. [I]. The two contributions from (F u — F d )Bi + 
(G u - G ,d )e 1 and from (Fg - F^)B 2 + (G u s - G d s )e 2 to T(B° d ) - T(B + ) are separately scheme- 
independent. Tab. U reveals that the former part is expected to give the dominant contribution 
to the desired width difference. The results also show a substantial improvement of the [X\- 
dependence in the NLO compared to LO. This dependence is plotted in Fig. [|for the two Wilson 
coefficients of the important vector operators. The approximation employed in pSp setting z = 
in the NLO correction is also plotted. Expectedly, the accuracy of this approximation decreases 
for small because the difference to the exact NLO result is of order a s (/ii) zlnz. For the 
final result of our coefficients we estimate the fii -dependence in a more conservative way: we 
vary fii in F u . . . Gg and F d . . . G d s independently. Further the variation with z and a s (M z ) in 
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z 




0.085 




0.070 


0.100 




m b /2 


m b 


2 m b 


m b 


m b 


^u,NLO _ pd,NLO 


0.865 
0.396 


0.270 
0.460 


-0.176 
0.386 


0.280 
0.469 


0.261 
0.452 


T-rtt,LO T-id,LO 
r s — r s 

t-,m,NLO T^d,NLO 


0.002 
0.035 


0.042 
0.033 


0.105 
0.026 


0.043 
0.031 


0.042 
0.035 


qu,NLO Qd,NLO 


-9.912 
-8.665 


-8.618 
-8.501 


-7.848 
-8.154 


-8.887 
-8.718 


-8.353 
-8.280 


s — s 

^u.NLO ^ycZ.NLO 


2.679 
1.668 


2.404 
1.850 


2.231 
1.902 


2.420 
1.854 


2.385 
1.843 



Table 1: Numerical values for the coefficients in @ for a s (Mz) = 0.118 and fx = m b = 
4.8 GeV. 



F u -F d 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

^/m b 



G u -G d 




-10.0 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

^/m b 



Figure 3: Dependence of F u — F d and G u — G d on [i\jm b for the input parameters in ( [R)| ) and 
fx = m b . The solid (short-dashed) line shows the NLO (LO) result. The long-dashed line shows 
the NLO result in the approximation of |]8p, i.e. z is set to zero in the NLO corrections. 



the ranges of (|16|) is calculated and all these sources of theoretical uncertainty are symmetrized 
individually and added in quadrature. The dependence on z is only an issue for G u — G d . We 
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find: 





NLO 


LO 


app 


jpu 




0.460 ±0.101 


0.270 ± 0.480 


0.440 ±0.119 


771 li 


-F* 


0.033 ± 0.046 


0.042 ±0.052 


0.025 ±0.045 


G u 


-G d 


-8.50 ±0.40 


-8.62 ±0.90 


-8.00 ±0.32 


G u s 


-G d s 


1.85 ±0.08 


2.40 ±0.23 


1.80 ±0.10 



The quoted central values correspond to the choice /xi = m b and the central values in dUSp. 
The third column in < JT7| ) shows the result for the approximation of [|8|], setting z = in the 
NLO corrections. For /ii = m b this approximation reproduces the size of the NLO corrections 
to F u - F d and G u s - G d s to better than 15% . The small NLO correction to G u - G d is, 
however, overestimated. The NLO result for this coefficient, which is largest in magnitude, is 
better reproduced by the LO result than by the approximation of [|]] . 

The origin of the z In z terms, which are the main cause of the discrepancy between 

the first and third column in flT7|), can be traced back to diagram D n of Fig. ^. This diagram 
defines the scheme of the charm-quark mass. One can absorb the a s (//i) z In z terms into the LO 
by replacing z with z = rnl(n) /m^(/x), which implies the replacement 

F u,(l) ^ F u,(l) _ (18) 

4:71 OZ 

in the NLO corrections to F u and similarly in the other Wilson coefficients. Here 7^ = 8 
is the LO anomalous dimension of the quark mass. This procedure sums the terms of order 
a"(/ii) z ln n z with n = 0, 1, ... to all orders in perturbation theory. This can be seen by per- 
forming an OPE of the transition operator T which treats m c as a light mass scale: then increasing 
powers of m c correspond to AB = operators of increasing dimension and m c and m b enter 
the result at the same scale [i\ at which the OPE is performed. In every order of the perturbation 
series In z is split into ln(/xf/m^) contained in the Wilson coefficients and ln(m^///f) residing in 
the matrix elements. Since there are no dimension-8 operators with charm-quark fields contribut- 
ing to r(£?°) — T(B + ), no terms of order ml ln(m^/^) can occur. From our NLO results we can 
indeed verify that the procedure in < JT8| ) removes the a s (ni) zlnz terms, while e.g. terms of order 
a s (ni) z 2 \nz persist as expected, because there are dimension-10 operators with charm-quark 
fields of the type m c (bq)(qb)(cc). Using z = 0.055 rather than z = 0.085 in the coefficients 
tabulated in the third column of ( |T7| ) indeed removes the disturbing discrepancy with the NLO 
result for G u — G d . Also the central values of F u — F d and G$ — Gg move closer to the NLO 
result, while no significant improvement occurs for F$ — Fg. 

The width difference in @ involves the product F qT B, which is independent of the renor- 
malization scheme and scales. In order to compare the scheme dependent coefficients F q with 
the calculation in [g] for z = 0, we need to take into account that the coefficients in [jSp are 
defined for matrix elements in HQET rather than in full QCD. The matching relation connecting 
HQET and full-QCD matrix elements of the four operators O used in []8p has the form 

(0) QCD (m b ) = (l + ^Mcpj (d) HQET (m b ) , (19) 
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where the 4 x 4 matrix C^ s can be found in Eq. (36) of The renormalization scheme of oper- 
ator matrix elements in full QCD is identical in our paper and in [|J§] . The only further difference 



is that the operators O are linear combinations, O 
with 



S 



i 

3 



4 
9 

V o 







2 


3 U 

o -I 
-1 o 



SQ, of our basis Q 
\ 



(Q, g s , T, T 6 



-4 





2 
3 



(20) 



/ 



(This simple relation holds beyond tree level because the renormalization schemes are identical. 
The preservation of Fierz- symmetry by the choice of evanescent operators in ([TJ]) is important 
for this property.) It follows that our coefficients F are related to the corresponding coefficients 
at scale \i = m b through 



A + ^B in 



3 V 4vr J 



at s + ^(b t s-a t c? is s 

47T V 



(21) 



Here we have suppressed flavour labels q = u, d and the double indices ij = 11, 12, 22 refering 
to the AB = 1 coefficients CiCj (see <JTT|)) . Note that in the notation of |||| labels u, d are 
interchanged with respect to our convention and that the coefficients with label 12 are defined 
with a relative factor of two. Using ( j2l"1 ) we have verified that the results of [|]] obtained for z = 
are in agreement with ours in this limit. 



3 Phenomenology 

3.1 r(B+)/r(BS) 

One can directly use (§) to predict the desired lifetime ratio: 

t(B+ 



i = t(b+) [m 

2 



r m- 



0.0325 



W, 



cb\ 



0.04 



B 



4.8 GeV/ \200MeV 



x 



1.0 ±0.2) B x + (0.1±0.1)B 2 - (18.4 ±0.9) ex + (4.0 ± 0.2) e 2 



(22) 



Here t(B + ) = 1.653 ps has been used in the overall factor and the hadronic parameters B\ . . . e 2 
are normalized at /x = nib throughout this section. 

In [J3p it has been noticed that without a detailed study of the hadronic parameters one expects 
t(B + )/t(Bj) to deviate from 1 by up to ±20%. This feature originates from the large coefficient 
of e% and persists in our NLO prediction in d22|), because the NLO corrections to G u — G d are 



small. Confronting (22) with the recent measurements [14, 1 



r(B° d ) 



1.082 ±0.026 ±0.012 (BABAR) 
1.091 ± 0.023 ± 0.014 (BELLE) 



(23) 
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one expects |ei| to be significantly smaller than 1/N C = 1/3, i.e. nonfactorizable contributions 
appear to be small. This result is confirmed by the existing computations of the e/s in quenched 
lattice QCD [Q,[T^]. However, due to its large coefficient sophisticated non-perturbative methods 



are definitely necessary to compute e% sufficiently accurately. The other important term in ( j22| ) 
is the first one: the NLO enhancement of F u — F d in ( |T7| ) has altered the coefficient of B\ in 



(g2j) from 0.6 ± 1.0 in the LO to 1.0 ± 0.2. While from the LO result not even the sign of this 
contribution was known, the NLO result now clearly establishes a positive contribution of order 
3% to t(£> + )/t(£>°) from the term involving B\. 

The hadronic parameters have been computed in [§] using the same renormalization scheme 
as in the present paper. They read 



(Bi, B 2 , e\, e 2 ) 



1.10 ± 0.20, 0.79 ± 0.10, -0.02 ± 0.02, 0.03 ± 0.01). (24) 



Using \ V cb \ = 0.040 ± 0.0016 from a CLEO analysis of inclusive semileptonic B decays []T6p, 
the world average f B = (200 ± 30) MeV from lattice calculations Q and m b = 4.8 ± 0.1 GeV 
in (122]), we find 



-(B- 



r(B° d ) 



1.053 ±0.016 ±0.017, 



1.041 ±0.040 ±0.013, (25) 



J LO 



where the first error is due to the errors on the NLO coefficients as given in ( J22| ) and the hadronic 
parameters (E3I), an d the second error is the overall normalization uncertainty due to m b , \V cb \ and 



Jb in (p2|). The first error reduces to 0.008 in NLO and 0.038 in LO, if the errors on the hadronic 
parameters are neglected, demonstrating the substantial reduction of scale dependence at NLO in 
comparison with the LO. This result is gratifying as the strong scale dependence observed at LO 
had been a major motivation for a NLO analysis. This is also seen in Fig. |], where we show the 
lifetime ratio as a function of the renormalization scale jii. We should, however, emphasize that 



the result and error given in (g5|) do not include the effects of l/m b corrections and unquenching, 
which could well be on the order of 0.05. The NLO result slightly exceeds the central value of 



the LO result and improves the agreement with the experimental value in ([231) . 



3.2 t(H°)/t(H 6 -) 

The SU(3)f anti-triplet (A;, ~ bud, S° ~ bus, ~ bds) comprises the ^-flavoured baryons 
whose light degrees of freedom are in a + state. These baryons decay weakly. Baryon lifetimes 
have attracted a lot of theoretical attention: the measured lifetime falls short of r(B®) by 
roughly 20% [^8j], which has raised concerns about the applicability of the HQE to baryons. Un- 



fortunately this interesting topic cannot yet be addressed at the NLO level, because t(A&) / t(B®) 
receives contributions from the SU(3)p-singlet portion T sing of the transition operator in (§) and 
NLO corrections to T sing are unknown at present. Further the hadronic matrix elements entering 
r (A;,) / t(B®) involve penguin contractions of the operators in (Q), which are difficult to compute. 
It is, however, possible to predict the lifetime splitting within the iso-doublet (H°, S^) with NLO 
precision. The corresponding LO diagrams are shown in Fig. [5|. For H&'s the weak decay of 
the valence s-quark could be relevant: the decays — > A b 7r~, — > A b e~V e and S° — > A fo 7r° 
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Figure 4: Dependence of t(B + )/t(B®) — 1 on /^i/m& for the central values of the input param- 
eters and no = m b . The solid (short-dashed) line shows the NLO (LO) result. The long-dashed 
line shows the NLO result in the approximation of i.e. z is set to zero in the NLO corrections. 



are triggered by s —> u transitions and could affect the total rates at the (9(1%) level [|19|]. Once 
the lifetime measurements reach this accuracy, one should correct for this effect. To this end we 
define 

r(5 6 ) = r(H 6 ) - T(~ b - A b X) = 1 ~ AbX) = -L- for ~ b = ~° b ,~ b , (26) 

where B(E b — > A b X) is the branching ratio of the above-mentioned decay modes. Thus T(E b ) 
is the contribution from b — > c transitions to the total decay rate. In analogy to (Q) one finds 

r(H 6 -)-r(H°) = G2 ?<\y 2 f B M B (\V ud \ 2 F u + \V cd \ 2 F c - F d ) ■ B S K (27) 

Here B~ b = (Li 6 (/x ), L^(/i ), L^O^o), L^K/io)) 7 comprises the hadronic parameters defined 
as 

(E° b \(Q--Q d )(fi )\E° b ) = f B M B M Eb Lf b (^), 
(E b \(Q u s -Q d s )(fio)\E° b ) = f 2 B M B M Eb Lf^ ), 
(E° b \(T u -T d )(^)\~° b ) = f B M B M Rb Lf{^ 
<H°|(T£-T|)(// )|H°> = f B M B M Sb L? s (v ). (28) 
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ws 



Figure 5: Weak scattering (WS) and PI diagrams for E b baryons in the leading order of QCD. 
They contribute to r(H°) and T(S^), respectively. CKM-suppressed contributions are not shown. 



In contrast to the B meson system, the four matrix elements in (|2g) are not independent at the 
considered order in kQcoj^b- Since the light degrees of freedom are in a spin-0 state, the matrix 
elements (E b \2Q q s + Q q \r> b ) and (H&|2T| + T q \E b ) are power-suppressed compared to those in 
(see e.g. [[]J(3J]). This, however, is not true in all renormalization schemes, in the MS scheme 
used by us 2Q q s + Q q and 2Tf + T q receive short-distance corrections, because hard gluons can 
resolve the heavy 6-quark mass. This feature is discussed in [0]. These short-distance corrections 
are calculated from the diagrams Ex . . . E± in Fig. |2|. For our scheme we find 



( L 1S 


[m b ) ^ 


= 


V L 2S 







a s [m b ) 

4:71 



-28/3 -7 
-14/9 7/2 



Lr b (m b ) 
L% b (m b ) 



(±QCD\ mC29) 



As an important check we find that the dependence on the infrared regulator drops out in (E^). 



With ( |29[ ) we can express the width difference in ( J27| ) in terms of just the two hadronic parameters 
Lf 6 and h\ h . We find 



- 1 



0.59 



W, 



cb\ 



0.04 



B 



T S 



-01 



X 



4.8 GeV/ V200MeV; 1.5 ps 
(0.04 ± 0.01) Lx - (1.00 ± 0.04) L 2 



(30) 



with Li = Lf b (fx = m b ). For the baryon case there is no reason to expect the color-octed matrix 
element to be much smaller than the color- singlet ones, so that the term with L 2 will dominate 
the result. The hadronic parameters L 1:2 have been analysed in an exploratory study of lattice 
HQET pD] ] for A b baryons. Up to SU(3)p corrections, which are irrelevant in view of the other 
uncertainties, Lf b and Lf b are equal. 



4 Conclusions 



We have computed the Wilson coefficients in the heavy quark expansion to order (AQco/mb) 3 
for the B + -B® lifetime difference at next-to-leading order in perturbative QCD. These coeffi- 
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cients depend on the scheme and scale ji^ used to define the matrix elements of the AB = 
operators in the effective theory. Our scheme is specified by the NDR prescription for 75, MS 
subtraction and the definition of evanescent operators given in ([T4]). The 0(a s ) accuracy is cru- 
cial for a satisfactory matching of the Wilson coefficients to the matrix elements determined with 
lattice QCD. Current lattice calculations, which are still in a relatively early stage in this case, 
yield, when combined with our calculations, t(B + ) /t(B%) = 1.053 ± 0.016 ± 0.017 [see ©]. 
The effects of unquenching and l/m& corrections are not yet included, but could well be on the 
order of 0.05. Next-to-leading order corrections to t(B + )/t(B®) were recently computed in the 
approximation m c = [§]. Taking the limit m c — > of our results we find agreement with this 
calculation. 

A substantial improvement of the NLO calculation is the large reduction of perturbative un- 
certainty reflected in the scale dependence of AB = 1 Wilson coefficients from the standard 
weak Hamiltonian. This scale dependence had been found to be very large at leading order, pre- 
venting even an unambiguous prediction of the sign of t(B + )/t(B®) — 1 up to now [Bp. With 
this major source of uncertainty removed by the NLO calculation, further progress will depend 
on continuing advances in the evaluation of the nonperturbative hadronic matrix elements and 
the computation of l/m b -suppressed effects. 



References 

[1] M. A. Shifman and M. B. Voloshin, in: Heavy Quarks ed. V. A. Khoze and M. A. Shifman, 
Sov. Phys. Usp. 26 (1983) 387; M. A. Shifman and M. B. Voloshin, Sov. J. Nucl. Phys. 
41 (1985) 120 [Yad. Fiz. 41 (1985) 187]; M. A. Shifman and M. B. Voloshin, Sov. Phys. 
JETP 64 (1986) 698 [Zh. Eksp. Teor. Fiz. 91 (1986) 1180]; I. I. Bigi, N. G. Uraltsev and 
A. I. Vainshtein, Phys. Lett. B 293 (1992) 430 [Erratum-ibid. B 297 (1992) 477]. For a 
recent review see: M. Voloshin, in: B physics at the Tevatron: Run-II and Beyond, Chapter 
8, [ hep-ph/020107Tl ]. 

[2] M. Beneke and G. Buchalla, Phys. Rev. D 53 (1996) 4991. 

[3] M. Neubert and C. T. Sachrajda, Nucl. Phys. B 483 (1997) 339. M. Beneke, G. Buchalla 
and I. Dunietz, Phys. Rev. D 54 (1996) 4419. 

[4] G. Altarelli, G. Curci, G. Martinelli and S. Petrarca, Nucl. Phys. B 187 (1981) 461. 

[5] A. J. Buras and P. H. Weisz, Nucl. Phys. B 333 (1990) 66. 

[6] Y. Y. Keum and U. Nierste, Phys. Rev. D 57 (1998) 4282. 

[7] M. Beneke, G. Buchalla, C. Greub, A. Lenz and U. Nierste, Phys. Lett. B 459 (1999) 631. 

[8] M. Ciuchini, E. Franco, V. Lubicz and F. Mescia, [ hep-ph/01 10375] ] . 

[9] D. Becirevic, [ |hep-ph/01 10121] ] . 



20 



The B + -B d Lifetime Difference Beyond Leading Logarithms 



[10] M. Di Pierro and C. T. Sachrajda [UKQCD Collaboration], Nucl. Phys. B 534 (1998) 373; 
Nucl. Phys. Proc. Suppl. 73 (1999) 384. 

[11] W. A. Bardeen, A. J. Buras, D. W. Duke and T. Muta, Phys. Rev. D 18 (1978) 3998. 

[12] S. Herrlich and U. Nierste, Nucl. Phys. B 455 (1995) 39. 

[13] M. Misiak and J. Urban, Phys. Lett. B 451 (1999) 161. 

[14] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 87 (2001) 201803. 

[15] K. Abe [Belle Collaboration], flhep-ex/0202009l ]. 

[16] D. Cassel [CLEO coll.], talk at Lepton Photon 01, 23-28 Jul 2001, Rome, Italy. 
[17] S.Ryan, [ pep-lat/01 1 lOlCQ . 

[18] D. E. Groom et al. [Particle Data Group Collaboration], Eur. Phys. J. C 15 (2000) 1; updated 
at http://pdg.lbl.gov. 

[19] M. B. Voloshin, Phys. Lett. B 476 (2000) 297. 

[20] M. Di Pierro, C. T. Sachrajda and C. Michael [UKQCD collaboration], Phys. Lett. B 468 
(1999) 143. 



